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Biquaternions algebra is complex expansion of quaternions algebra, which was elaborated 
by W.R.Hamilton in the middle of XIX century [1]. It is very useful and very suitable apparatus 
for description of many physical processes. In last decennial these algebras become be actively 
Ph ■ used in research of different problems of electrodynamics [2-4], quantum mechanics [5-9], fields 
C*~) ■ theory [10-13] and others. 

This algebra gives possibility to write many well-known systems of differential equations 
in theoretical physics in form of one biquaternionic equation in bigradiental form. Concept 
bigradient is theoretical generalization of gradient concept, differential operator which charac- 
terizes a direction and intensity of exchange of scalar field. By analogy bigradient characterizes 
a direction and intensity of exchange of biquaternionic fields which are used for description of 
complex scalar-vector physical fields. 

Here we elaborate the differential algebra of biquaternions on Minkowski space for construc- 
tion of generalize solutions of bigradiental biquaternionic equations, which are equivalent to 
the hyperbolic systems of differential equations of 8 order. We consider biquaternionic wave 
equation (biwave equation) and obtain on the base of Fourie transform its fundamental and 
generalized solutions. Shock waves are researched and conditions on their fronts are determined. 
£S) \ Invariance of biwave equation for the groups of orthogonal transforms, Lorentz and Poincare 
transforms are explored and relativistic formulas are constructed. The new scalar equations for 
potentials of bigradiental fields (KGFSh- equation)) is obtained, which contains Klein-Gordon- 
Fokk and Shrodinger operators, and its generalized solutions are defined. The biquaternionic 
representation of Maxwell and Dirac equations are considered, their fundamental and generalized 
solutions are obtained, which describe nonstationary, harmonic and static scalar-vector electro- 
magnetic fields, spinors and spinors fields. 

Properties of superposition of bigradiental operators, which is very easy calculated in biqua- 
ternions algebra, and methods of generalized functions theory make the process of solutions 
construction simple and fine, as you'll see hereinafter. 

1. Biquaternions algebra 

We introduce some notions and indications, which will be used here. Let e\, e^, e% are basis 
vectors of cartesian coordinate system in R 3 , eo = 1 ; F is three dimensional vector with 
complex components: F = F\e\ + F^e-i + F^e^ , Fj, f G C are complex numbers; Ejki is Leavy- 
Chivitta pseudotensor, 5jk is Kronecker symbol. 

We consider the linear space of hypercomplex numbers (biquaternions) B = {F = / + F} : 

aF + bG = a(f + F) + b(g + G) = (af + bg) + (aF + bG), Va, b G C, 

with operation of quaternionic multiplication [1]: 

FoG = (f + F)o( g + G) = (fg-(F,G)) + (fG + gF+[F,G}). (1) 
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Here and hereinafter we designate as (F, G) = FjGj the scalar product of vectors F and 
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G ; their vector product is [F, G] = Yl e juFjG k ei . 

j,k,i=i 

Biquaternions algebra is associative: 

FoGoH = (FoG)oH = Fo(GoH), (2) 

but not commutative: 

<F,G>=FoG-GoF = 2[F,G]. (3) 

It may be proved on the base of these properties of Leavy-Chivitta pseudotensor and basic 
elements: 

£jkl E-ljk £klj> E-jkl ^jlk E-kjh ^jkl^mnl 0~jmO~kn 0~jnO~kmi 

e o e = e , e oej = ej, ej o e k — -5 jk + e jk iei, 

(ej o e k ) oe m = -e jkm - 5 jk e m - 5 km ej + 5 mj e k = ej o (e k o e m ) , j, k, l,m,n = 1, 2, 3, 

( over similar indexes in product there are summation from 1 to 3, like to tensor convolution). 

From the properties of commutator (3) it follows 
Jacobi identity 

« F,G >,H > + « H,F >,G > + « G,H >,F >= 

= 4[ [F, G],H]+ 4[ [H, F],G}+ 4[ [G, H],F} = 

= -4(F(H, G) - G(H, F) + H(G, F) - F(G, H) + G(F, H) - H(F, G)) = 0. 

I.e. biquaternions algebra is Lie algebra. 

From (3) we see, that the product of two biquaternions is commutative if one of them is a 
scalar or if their vector parts are parallel. 

Some definitions. 

Bq. F - — f — F is called mutual for F = / + F . 

Bq. F = / + F is named complex- conjugate to F ( overline marks the complex conjugate 
number). IfFoF = FoF = l, then we name F unitary. 

F* = / — F is termed conjugate to F. If F* = F , then it is self conjugated. 
Selfconjugated biquaternion has the form: F = / + iF, where / and F are real. 

Scalar product of F 1; F 2 is the bilinear operation (F 1; F 2 ) = /1/2 + (-Pi, F 2 ) ■ 

Norm of F is nonnegative scalar value 



(F,F) = ^f.f+(F,F) = ^\f\'+\\F\\ 2 . (4) 
If F is real biquaternion (quaternion) then F* = F~ and 

||F|| 2 = F*oF = FoF*. (5) 

Pseudonorm of F is the value 



m) = Jf-f-(F,F) = J\f\ 2 -\\F\\ 2 , Re((F))>0. (6) 



It's easy to see, if F is self conjugated then F = F and 

FoF = FoF = ((F)) 2 . (7) 

If F o G = 1 then G is right inverse for F and it is termed as G = F -1 , thereafter F is 
left unverse for G and it is termed as _1 G. 
Simply the equalities are proved: 

(F o G)* = G* o F*, (FoG)" 1 = G" 1 oF" 1 . (8) 

Also simply the next lemma and theorem are proved. 

Lemma 1.1. // (F, F) ^ 0, then both inverse biquaternions exist and they are equal: 

F _1 = _1 F = F~/(F, F). (9) 

If (F, F) = then inverse Bq. does not exist. 

Theorem 1.1. By known F and B biquaternionic linear (bilinear) equations of a kind 

F o G = B and G o F = B 

have unique solution G = F" 1 o B and G = Bo F _1 appropriately, if and only if (F, F) ^ 0. 

Proof of this theorem follows from the lemma 1.1. 

If (F, F) = , then existence of the decision depends on B (see [14]). 

We introduce here the very useful Bq. 

3 = W + iP =^FoF* = (|/| 2 + \\F\\ 2 )/2 + i (lm(/F) + [Re F, ImF]) . 

((3)) 2 = W 2 - ||P|| 2 = HoE- 

Here W n P are real values. It describes the energy-impulse density of scalar-vector field in 
problems of theoretical physics. 

2. Lorentz and Poincare transformations on Minkowski space 

Let consider biquaternions on Minkowski space M = {(r, x) : r G R},x G R 3 } and groups 
of linear transformations on it. We quoternize M by using complex-conjugate Bqs.: 

Z = t + ix, Z = Z = t - ix, t G R 1 , x G R 3 . 

They are selfconjugated: Z = Z*, Z = Z* , and have equal norms and pseudonorms: 

||Z|| 2 = ||z|| 2 = r 2 + ||a:|| 2 = (Z,Z), ((Z)) 2 = ((Z)) 2 = r 2 - ||a:|| 2 = Z o Z (10) 

and 

Z- 1 = Z/((Z)) 2 , Z- 1 = Z/((Z)) 2 . 
Hence on light cone ( \t\ = \\x\\ ) inverse Bq. for Z, Z does not exist. 
Orthogonal transformation. Let consider conjugate Bqs.: 

U((p, e) = cos ip + e sin ip, U* = U _ = cos ip — e sin <p, \\e\\ = 1, ip G R 1 . 
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It's easy to see: 



(11) 



Lemma 2.1. Conjugate biquaternions U(y?, e), U*(</?, e) , tp E R 1 , define the group of 
orthogonal transformation on M, which are orthogonal on vector part of Bq. Z : 

z' = U o Z o u*, Z = U* o z' o u. 

This transformation is the rotation of space R 3 around vector e through angle 2(p . 

Proof. By calculating lemmas formulae we get conservation of scalar part and specified 
rotating of vector part: 

t' = t, x' — e(e, x) + (x — e(e, x))cos2ip + [e, x]sin2<p. 

Under (10)-(11) the pseudonorm ((Z')) 2 = UoZoTJ*oUoZoU* = ((Z)) 2 . Since r = t' , 
the norm of vector Z hold true: \\Z\\ = \\Z'\\ . 

Let consider superposition of two orthogonal transformation: 

Ui o U 2 = U 3 = u 3 + Us = (cos (pi + ei sin </?i) o (cos (f 2 + e 2 sin (p 2 ) = 
= cosy?i cos(/9 2 — (ei, e 2 ) sin^x siinp 2 + e 2 sirnp 2 cos ip 1 + t\ siny?! cos(/? 2 + 

+ [ei, e 2 ] sin y? 2 sin v? 2 . 

Since 

U3 1 = (U x o U 2 ) _1 = U 2 1 o u^ 1 = u 2 - o Uf, 

= (u x o u 2 )* = u;ou; = u 2 oUi= uj, 

we get: 

u 3 o u 3 = (Ui o u 2 ) o (u x o u 2 ) _ = Ui o u 2 o u 2 o ur = 1 

Hence U3 is also orthogonal transformation: 

U 3 = cosv? 3 + e 3 siny? 3! (p 3 = arc(cos u 3 ) , e 3 = ^ 3 

\\ U 3\\- 

Lemma has been proved. 

Lorentz transformations. Let consider mutual self conjugated Bgs. 

L(0,e) = che + ieshO, L = ch6 -ieshQ, 9 e i? 1 , ||e|| = 1 

(here hyperbolic sine and cosine are used). It's easy to see that they are unitary: 

LoL = ch6 2 - she 2 = 1. (12) 

Next lemmas are proved by use of simply calculations [14]. 

Lemma 2.2. Lorentz transformation has such biquaternionic representation: 

Z' = LoZoL, Z = LoZ'oL, ((Z')) 2 = ((Z)) 2 . (13) 
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Also easy to show that pseudonorm conserves on the associativity and unitarians (12): 

((Z')) 2 = LoZoLoL-oZoL- = ((Z)} 2 . 

If to enter designations: 

ch29 = (1 - v 2 y 1/2 , sh26 = v(l-v 2 y 1/2 , H<1, 

After calculating we get known formulas for scalar and vector parts of Z' and Z [13]: 

relativistic formulas 

/ r + v(e,x) (e,x)+vr 
t = — , x = [x — e(e, x)j + e — , 

V 1 — v 2 v 1 — f 2 

r = — , , x = (x' - e(e,x')) + e , y 

It corresponds to motion of coordinates system {Xi,X 2 ,X 3 } in direction of vector e with 
velocity v. 

Superposition of two Lorentz transformations with equal ej possess with group properties: 

Li o L 2 = L 3 = l 3 + L 3 = (chipi + ieshOi) o (ch9 2 + iesh9 2 ) = 
= {chQ x chQ 2 + sh6 1 sh6 2 ) + ie(sh9 2 ch6 1 + shQ x chQ 2 ) = 
= ch{6 l + d 2 ) + esh{6 1 + 6 2 ) 

But superposition of two Lorentz transformations with nonparallel e\ and e 2 does not constitute 
the Lorentz transformation: 

Li o L 2 = L 3 = l 3 + L 3 = (ch9i + ieishOi) o (ch0 2 + ie 2 sh6 2 ) = 
= {ch6\ch6 2 + (ei,e 2 )sh9ish9 2 ) + i(e 2 shQ 2 chQ\ + e\shQ\chQ 2 ) — [ei,e 2 ]sh9 2 sh9i 

What see, this Bq. contains real value in vector part. 
Poincare transformation. Note that 

L(0,e) = U(-i0,e). 

Using superpositions of these two transformations we get general form of the transformation, 
which may be named Poincare transformation. 

Definition. Poincare transformation on M is the linear transformation of the type: 

Z'=PoZoP*, Z = P oZ oP* , 

P = U o L = cos((p + iQ) + esin((p + i9), P* = L* o U* = cos((p — iO) — esin((p — id), 

which conserves the pseudonorm: ((Z)} = ((Z')). 

It's easy to prove because P o P~ = P* o P*" = 1 . By such transformations light cone 
(r = ||x||) is invariant set as ((Z)) = ((Z~)} = 0. 

Superposition of two Poincare transformations with equal ej also possess with group properties. 
It's easy to show: 

P 1 o P 2 = Li o Ui o L 2 o U2 = Li o L 2 o Ui o U2 = L 3 o U 3 = P3. 
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But superposition of two Poincare transformations with nonequal doesn't possess such 
properties. 

3. The space of generalized biquaternions B'(M) 

We will consider on M the functional space of Bqs. B(M) = {F = /(r, x) + F(r, x)} , where 
/ is complex function and F is three dimensional vector-function with complex components 
Fj , j — 1, 2, 3 . The partial derivative from Bq. over r or Xj we designate so: 

d T F = d T f + d T F, djF = ^- + — , j = 1,2,3. 

OXj OXj 

Let introduce two biquaternions spaces, basic one is 

B(M) = {$ = ^(r,x) + $(r,x)}, ^ G F>(R 4 ), G F>(R 4 ), j = 1,2,3, 

where D(R 4 ) is the space of finite infinitely differentiable functions on i? 4 , and conjugate 
space B'(M) = {F = / + F} of linear continues functionals on B(M) : 

3 

(F, *) = (/, ^ + ^(4^-), V$GB(M), 

i=i 

which be named the space of generalized biquaternions and we mark such Bq. with cap over it. 
Any regular Bq. F corresponds to the functional, which can be presented in integral form: 

(*..) = /(P(r,x), .(r,*» ^* 3 , V* 6 B(M). 

Bq. is singular if its action on B(M) can not be presented in such form. Example of singular 
Bq. is singular functions from D'(R 4 ) because D'(R 4 ) C B'(M) . Using them the more complex 
biquaternions may be constructed. 

In particular, in the mathematical physics problems the next singular functions are often 
used as simple layers. There generalization on B(M) are generalize Bqs. F5s , which define 
the functionals of the type: 

(F(J S ,*) = J (F(r,x),^(r,x))dS, V$ G B(M). 

s 

Here is the surface integral on a surface S C R 4 , which dimension may be equal 1,2,3. Defined 
and integrable on S Bq. F we name the density of simple layer as is customary. 

Differentiation. Analogically as in the theory of generalized function [15] using the definition 
of the partial derivatives of generalized Bq.: 

(djF, $) = -(F, dj&) for V$ G B(M) 

the derivatives of singular Bqs. and the derivatives of more higher order can be constructed. 

From the properties of differentiation of regular functions with finite gaps on some 3- 
dimensional surface S C R 4 there is 

a,-F = a,-F + n j [F] S (5s, j = r, 1,2, 3 (14) 
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Here djF is classic derivative, n = {n T , ni, n 2 , 113} is unit normal to S . In square brackets the 
gap of F on S stands, 

[F(t,x)] s = lim^ {F(r + en T , x + en) — F(r — en T , x — en)} , (r, x) E S. 

Definition. Bq. F is generalize solution of differential equation D(d T ,d x )F = G if 

(DF, = (G, $) for V$ G B(M) 

(here D is linear differential operator). If F is quasi everywhere differentiable regular function 
F then it is classic solution. 

Convolution of two biquaternions is the bilinear operation: 

3 

A(r, x) * B(r, x) = a * b - ^ (Aj * Bj) + (a* Aj) ej + (b * Bj) ej + eiji (Ai * Bj) e t , 

where usual convolutions of generalized functions [15] stand in brackets on the right. It's easy 
to see, here two operations of biquaternionic multiplication and convolution are united. 

In virtue of properties of convolution differentiation [15] we get the formula for derivative 
of Bqs. convolution: 

dj(A * B) = (djA) * B = A * djB, dj = d T , d u d 2 , d 3 . (15) 

It gives possibility to use more convenient form by its calculating . 

Fourier transformation. Generalized Fourier transformation (FTr.) of G is G which 
satisfies the equality: 

(G,*) = (G,*) for V$ G B(M). 

Here $ is classic FTr. of $ : 

mT , x)] = *(„,{) = J * (T , x)eMiT „ +i(i , x)) ir ixiix ^, * 6 B(M), 

which always exists in virtue of properties B(M) . If Bq. is regular and has classic FTr. then it 
is also generalized FTr. 

Also it's easy to prove (similarly in [15]), that 

F[A(r, x) * B(r, x)} = A(u, f) o B(u>, 0- (16) 

Using properties of FTr. of generalized functions we could get many similar properties of 
Fourier transformation for biquaternions, but we don't do this here. 

4. Bigradients. Poincare and Lorentz transformations 

Let consider special cases of differential operators which are representative for mathematical 
physics equations. But we will analyze them on B'(M) . 

We introduce operators, named mutual complex gradients: 

V+ = <9 r + iV, V" = 9 T -iV, (18) 



7 



where V = grad = (<9i, 82, 63) . We name them shortly bigradients. 

In the sense of given definitions their symbols are complex conjugate and selfconjugated: 
(V - )* = V~, (V + )* = V + . Their action is defined in according to biquaternions algebra 

V ± F = (8 T ± iV) o (/ + F) = (d T f =f i (V, F)) ± T F ± iVf ± i[V, F] 

(correspond to signs) or in conventional record 

V ± F = (d T f T i div F) ± d T F ± igrad f±irot F. 

It's easy to test that wave operator ( □) is presented in the form of superposition of mutual 
bigradients: 

□ = — - A = V~ o V + = V + o V", (19) 
or 1 

3 

where A = Yl djdj i s Laplace operator. Using this property, we can construct particular 
3=1 

solutions of differential biquaternionic equations on B'(M) of a type: 

V ± K(r,a;) = G(t,x), (20) 

which be named biwave equation. Solutions of this equations be named ± bipotentials of G . 

By Poincare transformations bigradients and biwave equations transform according consecutive 
affirmations. 

Theorem 4.1. By Poincare transformations (Z — > Z' = P0Z0P*) biwave equation 
transforms to biwave equation: 

D'K' = G', 

D = V ± = d T ± % grad x , D' = P _ oDoP = 9 r '±i grad x > , 
(relativistic formulas) 

K =P 0K0P, G =P 0G0P. 

Proof. See [14] and equalities 

D'K' = (P o D o P) o (P o K o P) = P o D o K o P = P o G o P = G'. 

Lets go to solving biwave equations. 

5. Generalized solutions of biwave equations. Shock waves 
Theorem 5.1. Generalized solution of biwave equation (20) can be presented in the form: 

K = V T G*^ + K . (21) 

where i/j(t,x) is simple layer on light cone r = ||x|| : 

^ = (47r||x||)- 1 5(r-||a:||), (22) 

which is the fundamental solution of wave equation: 

OJ, = 5(t)5(x). (23) 
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Bq. Ko(r,i) is a solution of homogeneous biwave equation (by G = 0) . It can be presented as 

K = V^{G *^o(r,x)}, (25) 
where ip (r,x) is the a solution of homogeneous wave equation: 

□V> = 0, (26) 

Mt,*) = J HO^P(H^^)±U\\r))dV(0, V0(O G MR 3 ), (27) 

Go(t,x) is arbitrary Bq., capable of this convolution with ip , or it can be presented as a sum 
of such type solutions. 

Proof. Using associativity and properties (15), (23) for first summand in (21) we get 
V ± K = V ± V T (G * V) = □ (G*V>) = G*O0 = G * 5(r)5(x) = G. 
Last equality in virtue of the property of 5 -function [15]. By analogy for the second summand: 
V ± K = V± V T {C *Mr,x)} = D{C * Vo(r, x)} = {□V'oM} * C = 0. 

Inversely, if K is the solution of homogeneous biwave equation, then every its components 
is the solution of Eq. (26). Thus and so it can be presented in the form (25) or decomposed in 
sum of forth such bipotentials for every components generated with different solutions of Eq. 
(26). In virtue of linearity their sum in (21) is a solution of Eq. (20). 

The form of solutions of wave equation (26), including fundamental one (22), is well known 
[15]. 

Fundamental solution of biwave equation we obtain if suppose G = S(t)S(x) . Then 
by use formula (21) we have its presentation: 

vI/(t, x) = d T ip ± i gradip. 
It gives possibility to construct solutions of biwave equation in form of the convolution: 

K = G * * + K 

for arbitrary right part of Eq.(20), accepting such convolution. 

Shock waves. Let consider generalized solutions of biwave equation. Whereas it is hyperbolic, 
there are solutions which are non differentiable on characteristic surfaces ( S ) where 

n 2 T = |M| 2 (28) 

(n T , rii, ri2, n$) is unit normal to S 1 in M, n — (ni, n2, n^), \\n\\ 2 = n\ + n\ + n§. This is a 
cone of characteristic normals of wave equation. 

Let define the conditions on the gaps of derivatives for regular solutions of biwave equation. 
By using formula (14) we get: 

V+K = V+K + {n r [k\ s -i (n, [K] s ) + n T [K] s + in [k\ s + i[n, [K] s }}5 s (t, x) = 
= V + K + (n T + in) o [K} s 5 s (t, x) = 
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Because V + K = then 

n T [k] s -i(n,[K} s ) = 0, n T [K] s + in[k\ s + i[n, [K] s ] = 0. (29) 

Such solution is named shock wave. In R 3 it has mobile wave front St with normal 
(711,712,713) , which propagates in virtue (28) with speed 

c=-tSt = 1- (30) 
||n|| 

This result is formulated in this theorem. 

Theorem 5.2. Shock waves satisfy to the conditions on its fronts : 

where m = n/||n|| is a unit wave vector directed along the speed of its expansion in R 3 . 
Proof. If to divide (29) per ||n|| subject to (30) , we get 

[k] 8 = -i(m,[K\ 8 ), [K\ 8 = im[k] 8 + i[m,[K\ 8 ] (31) 

Pass to its biquaternionic record we get the formula of the theorem. 
First equation (31) characterizes longitudinal shock waves. 

If to substitute it to second equation, we get formula for tangent component of vector K to 
wave front: 

[K] s - m ( m ,[K] s )=t[ m ,[K] s ]. (32)) 
They give the connection between gaps of real and imaginary parts of Bq. 

Generalized Kirchhoff formula for biwave equation Now we solve Cauchy problem 
for biwave equation. Let initial condition is given: 

K(0,x) = K°(x). (33) 

We have to construct the solution of Eq. (20) , satisfying to these data. 

We use for this generalized functions method. Let introduce regular functions G = H(t)G(t, x) , 
where H(t) is Heaviside function. If K° is regular then 

V ± K = G + 5{t)K°(x). (34) 

Therefore the solution is 

K(r, x) = V t {H(t)G * ip} + G(0, x) * ip + V T {K°(a;) * ijj}, (35) 

X X 

where sign * implies a convolution only over x . For regular G it has integral representation: 

X 

47rK(r,:r) = -V^j J G(r - r, y)^- + r' 1 J K°(y)dS(y)\ - r' 1 J G(0,y)dS(y). 

\r<r r=r ) r=r 

(36) 
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Here and further r = \\y — x\\, dV{y) = dy\dy 2 dy^ , dS(y) is differential of spheres surface 
r = r in point y . This formula is generalization of well known Kirchhoff formula for solution 
of wave equation [15]. 

6. Biquaternionic form of Maxwell equations and its modification 

As example of applications of biwave equation theory let consider Maxwell equation for 
electromagnetic field. It contains 8 differential Eqs. (2 vectorial and 2 scalar Eqs.). 
In space of biquaternions it has the form of biwave equation [2,4,16]: 

V + A + = 0. (37) 

Biquaternions of intensity A and charge - current of EM-field are 

A = + A = y/eE(T,x) +iy/JiH(T,x), = ip(r, x) + J(r, x), (38) 

where E is electric field intensity, H is magnetic field intensity; a density p and pulse J 
are denominated with use of density of electric charge p E and electric current density j E by 
formula: 

p = p E I yfe = y/edivE, J = y/pj E , 

e, p are the constants of electric conduction and magnetic conductivity; r = ct , t is time, 
c = 1/ y/epZ is light speed. 

Scalar and vector part of Eq. (37) gives us Hamilton form of Maxwell equations [17,18]: 

p = div A, d T A + % rot A + J = 0. (39) 

Real and imaginary part of Eqs. (20) gives us the system of 8 Eqs. of classical form. 

If to take mutual bigradient from Eq.(37), we get charge conservation law in scalar part and 
wave equation for intensities of EM-field in vector part: 

d r p + div J = 0, DA = irot J — gradp — d T J. 

Bq. energy-pulse of A-field 

S = W + iP=-AoA* (40) 
2 

contains energy density of EM-field 

W = 0,5\\A\\ 2 = ±{e\\Ef + p\\H\\ 2 ), 

and Pointing vector 

P = - [A, A] = c- l E x H. 
Relativistic formulas for Maxwell equations for Poincare transformation have the form: 
V+A' = 0', where A' = P o A o P, = P o o P. 

Shock EM-waves. From theorem 5.2 we get conditions on fronts of shock EM-waves: 

[A] s = im o[A] s . (41) 
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Its scalar and vector parts be of the form : 

(m,L4] 5 ) = 0, [A} s = i[m,[A} s ]. 
First condition shows that shock EM-waves are transversal. From second condition we have 

Vi[E] s = ^Jl[[H] s ,m], Jp [H] s = y/i[m, [E] s \. 
Using vectors of electric displacement and magnetic induction 

D = eE, B = pH, 

we state result. 

Theorem 6.1. On the fronts of shock EM-waves the gaps of intensities of EM- field satisfy to 
conditions: [E] s = c[[B] s ,m], [H] s = c [m, [D] s ]. That is equivalent to conditions: [D] s = 
c-^H^m], [B] s = c- 1 [ m ,[E] s ]. " 

Generalized Kirchhoff formula for Maxwell equation be of the form: 



47rA 



| J 0(r - r, y)^- + J A°(y)dS(y) I + r' 1 J 0(0, y)dS(y). (42) 

\r<r r=r ) r=r 



Integral representation for E,H can be get from here but they will be composite form. 

Modified Maxwell equations. In Hamilton form of Maxwell Eqs. (34) vector Eq. defines 
currents, scalar Eq. is definition of a charge. In bigradient of intensity of EM-field the scalar 
part is equal to 0. Corollary of it is law of charge conservation. I.e. Maxwell equations describes 
closed systems of electric charges and currents and generated by them EM-field. 

For open system the Maxwell Eqs. must be modified with use of scalar field a(r,x) in Bq. 
A = a{r,x) + A(t, x) . Then from Eq. (32) follows its modified form: 

d T A + i rot A + J = grad a, p = div A — d T a. (43) 

If p and J are known this system for determination a and A is closed. Scalar field a(r, x) 
may be called resistance- attraction field [11,12]. 

Comment. Biquaternionic form of Maxwell equation (32) has deep physical sense: electric 
charges and currents are physical appearance of EM-field bigradient, and intensities of EM-field 
are the bipotential for its charge-current. 

7. Biquaternionic presentation of Dirac operators and their properties 

Biwave equation (20) cam be written in the matrix form: 

3 

J^D^jbj^gm, m,j = 0,1,2,3, (44) 

3=0 

where 6 = b, go = g, bj = Bj, gj = Gj, j = 1,2,3; and D^j are the components of matrix 
D ± (corresponding to sign), which have the form: 

d T —id\ —id 2 —ids 1 ( d T id\ id 2 id 3 

ir = D= { Wl ^ -f ld l , D- = D= \ % f ~f ) (15) 

S ld 2 10 3 O r —lO\ I I —IO2 — 2C>3 O t lOl " 

id 3 —id 2 id\ d T J [ — id 3 id 2 —id\ d T 
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Easy to test that their product (operators superposition) satisfy to relation: 
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3 

^2 D mj D ji = SmiH, j,m,l = 0,1,2,3. (46) 

3=0 

We'll show that matrixes (45) are the differential matrix operators of Dirac, which possess by 
such properties [13,17]. For this we afford them in matrix form: 

3 

D + = J2^d 3 , (47) 

3=0 

where, as follow from (45), matrixes have the components: 

-i ] ( -i 

D o = I Dl = < i 2 = i . 3 

/ [' i % < 

-% J [ -i 

here / is unit matrix. As you see, this is 4-dimensional Dirac matrixes, composed from 2- 
dimensional Pauli matrix: 

1 \ ( ±i \ ( ±i \ f ±i 
OlJ'VO Ti J ' \±i J ' \Ti 

Therefore Dirac matrix D ± has biquaternionic presentation V ± . 
The differential operators: 

D+=V + + m, D-=V"+m 

constitute biquaternionic presentation of matrix Dirac operators: D ± + ml . By m = it is 
biquaternionic form of Maxwell operator. By this cause biquaternionic equation of form: 

D^B = (V ± + m) o B = F, m e C, (48) 

is the biquaternionic presentation of generalized Maxwell-Dirac equation (MD- equation). 
After simply calculation we get this lemma. 
Lemma 7.1. 

D+oD" =D" oD+ =n + m 2 + 2md T . 

For imaginary m = ip 

r>i p or>; p = n-p 2 + 2i P d T . (49) 

It is interesting that right part (49) contains Klein- Gordon- Fokk operator (□ — p 2 ), and 
Srodinger operator ( A + 2ipd) T . By this cause we will name the equation 

Du + 2md T u + m 2 u = /(r, x) (50) 

Klein- Gordon-Fokk-Srodinger equation or shortly KGFSh-equation. 

As you'll see further, addition of 2md T to KGF-operator essentially simplify construction 
of solutions of Maxwell and Dirac equations and their representation. 
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8. Generalized solutions of Maxwell-Dirac equation (48). 
KGFSh-equation and scalar potentials 

Using lemma 7.1 and theory of generalized functions it is easy to construct generalized 
solutions of MD-equation(48). 

Theorem 8.1. Generalized solutions of MD- equation (48) have the form: 

B = B° + o (ip*F) = B° + * * F, (51) 

where B°(r, x) is a solution of homogeneous equation: 

D± o B° = 0, (52) 

is its fundamental solution: 

Vf^r, x) = d T ip + mip ± i gradip, 

ip(r,x) is fundamental solution of KGFSh-equation(50). 

Proof. In virtue of linearity of equation, it is sufficiently to prove this assertion for the 
second summand in formula (51). If to substitute it to Eq.(48), by use of properties (15) and 
5 -function, we get: 

D±D^ (ip* F) = F * (D0 + 2md T ip + m 2 ip) = F * S(r)S(x) = F. 

Ex facte, whatever solution of Eq. (48) can be presented in form (51). 
Let formulate this theorem for imaginary m — ip . 
Theorem 8.2. Generalized solutions of MD-equation : 

(V ± + ip) o B = F (53) 

can be presented in the form: 

B = B + (V T + ip) o (F*V>) , (54) 

where ip(r,x) is fundamental solution of KGFSh-equation: 

D^-p 2 ^ + 2ipd T ij = 5(r)5(x), (55) 

B (r, x) is a solution of Dirac equation: 

(V ± + ip) o B = 0. (56) 

From these theorems it follows that solutions of KGFSh-equation determine solutions of 
MD-equation. 

Theorem 8.3 Generalized solutions of KGFSh-equation (50) are presented in the form: 
u = f * ip + uq , where ip is its fundamental solutions: 

iP = ( ae - m|WI 5(r - \\x\\) + (1 - a)5(r + ||rr||)e m||a;|1 ) + rfo, Va G C, 

t:7T uXM 
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5(r± \\x\\) is simple fiber on light cone \\x\\ = =)=r , uo(t,x) is a solution of homogeneous Eq., 
which exists only for imaginary m = ip and 

«o(r,x) = e-^ J ^(Oexp(t(($,x)±||e||r))dK(0, V0(O G Li(i? 3 ). (57) 

Proof. From equation for fundamental solution: 

+ m 2 ^ + 2m<9 r ^ = S(t)S(x), (58) 

follows Helmholtz equation: 

{A-k 2 }F T [ip] + 5(x) = 0, k = iu-m. 
fundamental solutions of which are well known: 

FrbP] = T-ir^r (ae {tuJ - m)llxl1 + (1 - a) e - ( ^- m)INI ) , Va G C. 

47T ||x|| 7 

From here by use of properties of direct and inverse Fourier transformation the first formula of 
theorem follows. 

Support of first summand is expanding over time sphere of radius r ( t > ) , and for 
second one it is converging over time sphere of radius |r| (t < 0). 

If m is imaginary ( m = ip ) and support of solution over time is r > , then 

* = 5rw ,(T " W) - (59) 

It's interesting that here density of simple fiber in light cone is fundamental solutions of 
Helmholtz Eq. with wave number p . 

Solutions of homogeneous KG FSh- equation. Let consider 

Ou + m 2 u + 2md T u = (60) 
which in Fourier transform space satisfy to algebraic Eq. 

(ll£llM^ + ^)>>,0 = o. (6i) 

If Rem ^ 0, then ||£|| 2 — (u + im) 2 ^ by V{£, u} G R 4 . In this case the equation have only 
trivial solution: u* = . 

Provided m = ip this Eq. has infinitely set of singular solutions of the type: 

«>,0 = ^,0«(liei| 2 -(w-p) 2 ), (62) 

where (f>(u,£) is a density of simple fiber. It is arbitrary given and integrable function on cones 
||£|| = |w — p\ - By calculating of originals, we get: 

oo 

u(t,x) i= i J du J exp (-«(£,£) -iur)dS(£) = 
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{4>(f> +UL «e"" lfl|T - 4>(f> - K\\ , i y m,T } exp x)) rf6<*6<i&. 



In virtue of arbitrary 0, we get formula (57). The theorem is proved. 

Harmonic vibrations. Let construct solutions of MD-equation (48) in case of harmonic 
vibration with frequency uj : 

B = B{x)e iuJT , F = F{x)e iulT . 
Then we have the equation for complex amplitude of type : 

V± +p) B(x) = (u + p + V) o B(x) = F(x). (61) 
which be named gradiental. It easy to calculate 

V (w) oV(_ w ) = -(cd 2 + A). 
By using this property it easy to prove the theorem. 

Theorem 8.4. Generalized solutions of gradiental equation (61) can be formulary: 

B = Vj +p (x* F) + Sp^\ 

where 

ae -*fc|MI ( a _ i) e <*l|x|| 

x = -17g + 4, Ml ■ * = l- + "l^°- VaeC - < 62 > 

Sp(» + "> = VJ +p (xo*C(x)), 
Xo(x) a solution of Helmholtz equation with wave number k: 



Xo (x)= J p(e)e- ik ^dS(e), 



e =1 



p(e) is V integrable on unit sphere function, C(x) is V Bq., for which this convolution exists. 

If u > first summand in (62) describes radiated spherical wave, second one - converging 
spherical wave. The wave number k defines the length of these waves A = 2n/k, which (for 
positive p) decreases by increasing p. But if p < then A in interval {— w,0} increases by 
increasing \p\ and by p = —u spherical waves vanished. 

Statics. This theorem also defines solutions of static MD-equation (when u — ) . 

Theorem 8.5. Generalized solutions of static MD-equation can be formulary: 

B = VJ ( X *F) + Spto 

where 

ae'^W (a - l)^" 

47T ||x|| 47T ||x|| 

S P (P) = VJ (xo*C(x)), 
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Xo(x) is a solution of of Helmholtz equation with wave number \p\ : 

Xo(x)= J p(e)e- i ^dS(e), 

||e||=l 

p(e) is V integrable on unit sphere function, C(x) is V Bq., for which this convolution exists. 

We remark that Bq. C G B'(M) . For regular B this formulas gives classic solutions of 
MD-equations. 

If p = from these two theorems we get solutions of modified Maxwell equations in cases 
of harmonic vibrations and statics. 

9. Spinors and spinors fields 

In quantum mechanics solutions of Dirac equation: 

(V ± + ip) o Sp = 0, Rep = 0. (63) 

are called spinors [13,17,19]. 

Theorem 9.1. Spinors may be formulary: 

Sp = D5o(^*C(r,x)), (64) 

where ipo is a solution of homogeneous KGFSh-equation, or can be presented as the sum of 
such form solutions, VC G B'(M) for which this convolution exists. 

Proof. By substituting formula (64) into Eq.(63) we get: 

D± Sp = D±D^ (fo * C) = (Ctyo + 2mc^ + m 2 ^) * C = 0, m = ip. 

Inversely, if spinor (64) is a solution of (63), then 

(□ + 2md T + m 2 ) o Sp = D^D± Sp = D^ t = 0. 

Therefore scalar part and components of vector part Sp are the solutions of KGFSh-equation 
and Sp may be formulary as the sum of solutions of type (64). 

In particular case, when C = 1 , we have spinor of scalar ip -field : 

*o = ( V=F + m) o^ = mipo + d T i>Q =F % grad^ , 

Consequently C-field spinors , generated by potential ip Q , may by formulary: 

Sp ± = *±*C(r,a;). (65) 

Scalar harmonic potentials. Let see formula (57) where there are two plane harmonic 

wave: 

(pf(r, x) = exp (i ((£, x) - pr ± ||f || r)) . 

They also are solutions of homogeneous KGFSh-equation. Wave vector £ defines direction 
of wave motion. The length of these waves is equal to A = 2tt / ||£|| , their frequencies are 
u> — \p ± ||£||| , periods T = 27r/ \p ± ||£||| . In depend of sign, one of them has the supersonic 
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phase speed ( V > 1 ) , other one has the subsonic phase speed ( V < 1 ) , as V — 1 ± . 
V — > 1 ± by ||£|| — >■ oo , but ui — > oo . By ||£|| — > \p\ velocities V — > 1;0 and frequencies 
u) — > j } ; oo (corresponding to sign). 
Generated by these waves spinors 

{VT+ip)o^{T,x)=±{lU\\+i)^- 

have the form, norm and pseudonorm: 



S = Sp± o (Sp±)* = 1 - =► ||E|| = 2, ((H)) = 0. 



||Sp±|| = l, <(Sp±» = 0. 
We name them elementary £ - oriented harmonic spinors. Its energy-impulse is equal to 

n 

C -field £ -oriented harmonic spinor is the spinor: 

Cp 5 = C(t,x) *Spf(r,x). (66) 
From formula (64) it follows that C -field harmonic spinors may be formulary as 

Cp = C(t,x) * J Sp±(T,:r)0(£Ki«3, 0(0 e ^i(^ 3 ), 

Elementary harmonic (cu + p) -spinors are defined as: 

^ +P) (x, e) = ( V + to + p) o e -^ e ' x > = _L ( w + p - ifce) e"**^) (68) 
A;V2 

It easy to test that 

||^ + "|| = 1, = 0. (69) 

Here e is the direction of the spinor, k = \to + p\ is its wave number. Energy-impulse of ^q +p) 
are equal to 

S = ^ +p) o {^ +P) Y = 1 - i e stgn(co + p) 
Theorem 9.1. G-field harmonic spinors can be formulary 

Gp w+P (x, e) = G(x) * &q +p (x, e) (e-oriented spinors field) 

" T a^,,-Q(„./^(«.. W .) (— mm 

!|e||=l 

Vp(e) G Li({e G i? 3 : ||e|| = 1}). 

Static spinors are obtained by u> — . 
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Conclusion. The solutions of biquaternionic form of Maxwell-Dirac equations are here 
received in class of generalized function. That allows to build the decisions as for regular 
biquaternionic functions, so and at presence of singular sources in its right part. It's possible to 
use at building of biquaternionic theories of the elementary particles. At calculation of spinors 
fields, it is possible to flip differentiation on components C-field, when this suitable. C-field 
too can be singular generalized function. 

Bigradients, biwave equations and their decisions were used by author earlier for building 
of one models of electro-gravymagnetic fields and interactions [11,12,20]. It's possible to find 
much other useful applications of the differential algebra of biquaternions, what is offered to 
interested reader. 

Key words: biquaternion, bigradient, biwave equation, Maxwell-Dirac equations, generalized 
solutions, shock waves, spinors, harmonic spinors, spinors field 
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